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ABSTRACT 

Equations governing the motion of a specific class of singularities of the Euler equation in 
the extended complex spatial domain are derived. Under some assumptions, it is shown how 
this motion is dictated by the smooth part of the complex velocity at a singular point in the 
unphysical domain. These results are used to relate the motion of complex singularities to 
the stability of steady solutions of the Euler equation. A sufficient condition for instability 
is conjectured. Several examples are presented to demonstrate the efficacy of this sufficient 
condition which include the class of elliptical flows and the Kelvin-Stuart Cat’s Eye. 
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1. INTRODUCTION 


Singularities of the Euler equation for inviscid incompressible flow have been studied ex- 
tensively for many years. An important question that remains quite controversial is whether, 
in three dimensions, an Euler flow can develop a singularity in finite time when the initial 
conditions are smooth (see Pumir & Siggia 1 , Majda 2 , Kerr 3 , Brachet et al 4 and the references 
therein). While this paper does not address this controversy directly, except within a narrow 
class of singularities under restrictive assumptions, we expect that the approach presented 
here may help in a fuller resolution of this problem. All previous singularity studies have 
been limited to examining solutions to the Euler equation in the real (physical) spatial do- 
main, though in some cases the singularity presence in the complex unphysical plane was 
deduced from monitoring the Fourier coefficients. 

For simpler Euler flows with sharp discontinuities like the Kelvin-Helmholtz problem, 
previous work by Moore 5,6 , Caflisch & Orlenna 7 and Cowley et al 8 have shown the advan- 
tage of studying the equations in the complex unphysic.al plane. In that case, the singularity 
formation in the physical domain can be understood in terms of two processes: (a) instanta- 
neous singularity formation at some points in the complex unphysical domain at the initial 
time, and (b) the motion of such complex singularities towards the physical domain. The 
question of finite time singularities in the real domain then becomes a question of whether 
such singularities approaching the real domain actually impinge it in finite time. Moore 5 ' 6 
and Caflisch & Orlenna 7 have shown that for specific symmetric initial conditions, the sin- 
gularities form on the imaginary axis after a finite time, which then move towards the origin 
leading to a singularity formation on the real axis. However, recent work by Cowley et al 8 
suggests that the Moore singularity 5,6 actually forms at the initial instant of time at points 
off the imaginary axis. Later on, some of these singularities merge in pairs on the imaginary 
axis before traveling towards the origin (which is part of the real domain). Indeed, it appears 
that in a wide class of differential equations including the familiar one-dimensional inviscid 
Burger’s equation, singularities that show up later in the real domain are either present at 
some complex location initially or are formed at the initial instant of time. Such singulari- 
ties typically form at points (including infinity) in the complex spatial plane where a regular 
perturbation series in time (for early times) becomes disordered. In some cases, it is found 8,9 
that singularities are born at infinity in the sense that the singularity location in the finite 
spatial plane recedes to infinity as t — > 0 + . The study of complex singularities appears to 
have been a fruitful exercise that transcends the question of finite time singularity. Indeed, 
recent work on an interfacial evolution problem by Tanveer 10 suggests that this may be an 
important step in understanding the effect of small regularization on an otherwise ill-posed 
time-evolving flow in the physical domain. 
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Despite successes in each of the cases cited above, and perhaps in other cases unknown 
to the authors, the problems were reducible to one independent space-like variable besides 
time. We are unaware of attempts to generalize these ideas for multiple spatial variables 
as appears necessary for studying a general Euler flow. Here, we study a class of moving 
singularities of the Euler equation whose form appears to remain unchanged with time, with 
the expectation that, for certain other initial conditions, the process (a) described in the last 
paragraph leads to the formation of complex singularities of this form when they are not 
present initially. 

In this paper, we provide analytical evidence to show how each singularity of a specific 
type moves with the smooth part of the complex velocity field evaluated at that point. Of 
course, in a simpler vein, it has long been known that point vortex singularities in the physical 
domain move with the local velocity at the singularity location 11 . However, with certain 
assumptions, this result is generalizable to singularity motion in the complex unphysical 
plane as well. This result is used to conjecture sufficient conditions for the instabdity of a 
steady solution of the Euler equation. We present two conjectures. Conjecture 1 addresses 
the absolute instability of any steady Euler flow with a stagnation point. Conjecture 2 
addresses instability in the Langrangian sense; when the conditions of this conjecture hold, 
there exist disturbances that grow as we follow some fluid particle. 

Several applications to the hydrodynamic stability of inviscid flows in an infinite or peri- 
odic domain — which include the class of elliptical flows and the Kelvin-Stuait Cat s Eye 
are made to illustrate the use of Conjecture 1 to determine the absolute instability of these 
flows. 

2. MATHEMATICAL EQUATIONS IN THE COMPLEXIFIED DOMAIN 


Consider the Euler equation for the motion of an invisicid fluid where the vorticity u>(x, t ) 
and the velocity u(x,t) are determined by 



u> t + u ■ Vcu — • Vu = 0 

(1) 


V • u = 0 

(2) j 

■ 

V xu = w 

(3) | 

i 


Each of the components of the above equations can be put in the form 

| 


F(xi,x 2 ,x 3 ,t ) = 0 

(4) | 

s 


where x = (zi , x 2 , x 3 ). In the case of 27 t periodicity in each direction, if we assume that each j 

“I 

coefficient Cju in the following representation of F 



F(x„ l2 ,x 3 ) = ZC i u^ t ^ 

(5) ! 

i 


hk,l 
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satisfies the relation 


\C jU \ < K e~ s M +w+m (6) 

for some constant K and 8 which are independent of j , k and /, then it follow's that (4) will 
hold for complex x\, x 2 and x 3 when | Im xj|, j/m x 2 \ and | Im x 3 | are each small. Since the 
right hand side of (4) is zero, which is an analytic function of x\, x 2 and x 3 in the complex 
domain, it follows that (4) holds for all complex x. The condition (6) will certainly hold for 
the Fourier series coefficient of each component of equation (1), (2) and (3) if each component 
of u and u> and its derivatives satisfy such a condition. We restrict our attention to cases 
when this is true at some early time, which we assume is possible for a certain subclass of 
initial conditions. Then we arrive at the conclusion that each of equations (1), (2) and (3) 
hold for complex x (i.e. each component is complex ). However, since the right hand side 
of each component of (1), (2) and (3) is zero, which is analytic for all complex x and t, it 
follows that (1), (2) and (3) will hold for all x and t. The singularities of u and uj are those 

that are consistent with the solutions of (1), (2) and (3) in the complex x domain. The. 

same arguments hold even when the periods in xj , x 2 and x 3 are not the same or when the 
geometry is infinite in some directions and periodic in others. 

We examine a certain kind of initial conditions for which the velocity and vorticity are 
each real and analytic everywhere for real x, but have complex singularities at a singularity 
surface determined by 

d{x, 0) = 0 (7) 

where d is an analytic function of each component of x that is real valued and greater than 0 
for real x. Here d is also a function of t and its time evolution equation will be derived later. 
We also require the non-degeneracy condition that 'S/d ^ 0 where d = 0*. We consider 
initial conditions that can be decomposed into 

u(x, 0) = u s (x, 0) + /(d(x, 0)) q{x, 0) (8) 

u(x, 0) = w 5 (x, 0) + /'(d(x, 0)) p(x, 0) (9) 

where 

f(d) = d° (10) 

with 0 < a < 1. Thus the initial vorticity blows up at the complex singular points, though 
the velocity is finite. In (8) and (9), we assume that each of u s (x, 0), w 5 (x, 0), <f(^>0) and 
p(x,0) are analytic functions of x and that they are real valued for real x. Thus the initial 

*We have shown later that this condition is satisfied at later times once satisfied initially. 
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conditions (8) and (9) correspond to smooth initial conditions in the real domain. We also 
require that 


V • u s (x, 0) = 0 

(Hi 

t 

V x u s (x, 0) = cD s (x, 0) 

(12) 

( 

} 

and 


1 

V x [fq(x,0)] = f'p{x,0) 

(13) 

■ 

Now, for t > 0, we seek solutions to (1) - (3) of the form 


| 

] 

u(x,t) = u s (x,t ) + f(d(x,t )) q(x,t ) 

(14) 

■ 

^(S,t) = a s (x i t) + f\d(x,t))p(x,t) 

(15) 

m 

From (2) and (3), we get the following relations: 



V •«, + /' jjV-q + q-Vd =0 

(16) 


f 1 

V X u s + f —V X <f + Vd X <f — uj s + f'p 

(U) 


If we now require that 



V • u s = 0 

(18) 


and 



V x u s = C3 S 

(19) 

= 

u> St + u s • V LO s U) s • = 0 

(20) 


then the relations (16) and (17) imply that 


: 

L\7-q + q-Vd = 0 

(21) 


X q + Vd X q = P 

(22) 

= 

In addition, the substitution of (14) and (15) into (1), with the use of relations (20) 

■ (22), 


yields 


i 


f"p[d t + u s ■ Vd] + /' p t +u s • Vp — p^-p-V • q + f{q- V) p + - (<7 • V) 

■ V d) - (p • V) u, - i (u5 a • V) <?- ,7(V<1) '(Vxfl- /p • V,- = 0 (23) 
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Now notice that if we choose to evolve d(x, t ) according to 

dt + u s • V d = 0 (24) 

then 

Pt + Us- Vp- ■ q + /(<f • V) p + -^(<f • V) w s = <f(tu 5 • V (1?) 

+ (P ' ’ V) (u5 a • V) q + qf(Vd) -(Vxq) + fp • V<f (25) 

The existence of p satisfying (21), (22) and (25) for small enough time follows from the 
existence of u5 and u; 5 , each satisfying (1) - (3), and of d satisfying (24), since tu — to s = f f (d)p. 
We also notice that the coefficients in (21), (22) and (25) are continuous even at d = 0. 
Therefore, it appears reasonable to assume that for small enough time, with appropriately 
regular initial data, such a solution is C l in X\, X 2 and at d = 0 in the extended complex 
domain. Since, (24) implies that d > 0 in the real domain, once it is so initially, it follows 
that such a solution p will be real and analytic in the real domain for sufficiently small t 
since it is so initially. Notice that if we were to choose d t + u s • Vd / 0 at d = 0, then the 
resulting evolution equation for p would contain singular coefficients in x at the singularity 
surface d = 0, which would preclude a C 1 solution for p at d = 0. Indeed for a general choice, 
p is not even expected to be finite at d = 0. 

Despite the choice (24) for the evolution of d - which results in (25) having C x coefficients 
when d = 0 - we cannot rule out the possibility that the system of nonlinear evolution 
equations for cj 5 , p and d develops spontaneous singularities after a finite time. Hence, 
global existence in time cannot be expected in general. However, we assume that there 
exists a nontrivial subclass of initial conditions so that no spontaneous singularity develops 
in the complex domain for p, <f, u s or o5 s , other than a weak singularity of p and q at d = 0 
(in the sense that p is C x but not analytic). In that case, the continuity of p at d = 0 implies 
that at later times, the singularity of Q is still determined by 

d(x,t) = 0 (26) 

where f* is singular. Writing (24) in the characteristic form, it is clear that (26) will be 
satisfied for points Xq ( t) determined by 

^ = «,(x „((),«) (27) 

with the initial condition xo(0) satisfying the constraint that 

d(x o (0), 0) = 0 (28) 
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Thus, each point on the singularity surface moves with the speed u s evaluated at that point. 

However, u s is the smooth part of the vector field u at d = 0 in the complex x plane, which 
means the location of complex singularities of the type considered here can be determined 
without having to calculate /T(;r, i) from (25). Furthermore, with the assumption on the j 

smoothness of u s and continuity of p for all time, any singularity of this kind can only « 

approach the physical domain without actually impinging it in finite time. This can be 
deduced from (24) by putting it in the characteristic form and noting that u s is real in the 
real domain and d(x, 0) > 0 in the real physical domain. Furthermore, if a flow has a real 
point (say the origin, without loss of generality) where u s - 0 with a velocity gradient tensor 
T kj = independent of time, then replacing u s in (27) by T x shows that a complex 
singularity can only approach the physical domain exponentially in time. This observation = 

may be relevant to the numerical calculation of Hrachet et al in the Tailor-Green problem, : 

where they note an exponential approach of a singularity towards the physical domain. : 

However, since the deformation tensor is traceless, because of the divergence condition, it is | 

clear that only a restricted set of points on the initial singularity surface are a subset of the i 

complex stable manifold of the fixed points | 

Equation (27) is convenient since it shows that one need not solve for p(x, t) to determine I 

the singularity location in the complex plane even though one needs to determine u s (x } t) 

in the complex plane. When exact solutions are known for which may be steady or f 

unsteady, this poses no problems; however, when the solution is only determined numerically i 

in the real domain, one cannot expect to determine u s in the complex plane through Fourier j 

extrapolation, since such a procedure is ill-posed. In such a case, the solutions have to be i 

determined numerically, directly in the complex plane. This will be the subject of another | 

paper. However, if the singularities are close to the physical domain, it is possible to obtain I 

complex singularity trajectories directly from u s in the real domain. This can be seen by I 

decomposing xq (/) = #(/) T i?/(f), into its real and imaginary parts and taking the real and | 

imaginary part of (27) in the limit of small y. We then get I 

§ = 3,(i(<),<) ( 29 ) 

at 

^ = Vu s (i(t),t) y (30) 

at 

Equations (29) and (30) are the approximate singularity trajectory equations when a complex 

singularity is close to the physical domain. In this case, one needs ii s and its gradient in only - 

the real domain, which can be determined numerically. 

§ Complex stable manifold of a real fixed point of (27) is defined as the set of complex initial conditions 
x o (0) such that xo(/), satisfying (27), tends to the fixed point as t — * oo. 
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Once the singularity locations are determined, one can determine p(x 0 (t), t ) and Vd(xo(t ), t.) 
in terms of a relatively simple system of ordinary differential equations, as will be derived 
shortly. These quantities are of interest since for a point x sufficiently close to Xo(t) (so that 
d is small) but not on the singularity surface, we have 


J ~ a[Vd(xo{t) y t)-(x-x 0 {t))}- { '- a) p(x 0 (t),t) 

(31) 

For brevity of notation, we define 


A(t) = Vd(x 0 (f),f) 

(32) 

B(t) = q(x 0 (t),t) 

(33) 

C(t) = p(x 0 (t),t) 

(34) 

Then, it is clear that, since d(xo(t),t) — 0 and each of /, f /f are zero at d = 
and (22) implies that 

0, then (21) 

AB = 0 

(35) 

Ax B = C 

(36) 

It follows that 


C x A = (A- A)B 

(37) 

Furthermore, it is clear that (25) implies that 


f -(4-^ + xtJ 

(38) 

where T is a second rank tensor whose elements are defined by 


T - dUs > 
T>k - 

(39) 

Assuming A • A / 0^ we get from (37) that 


dt 

(40) 

where S is a second rank tensor whose elements are given by 


— * 

w s • A 

Sjk — Tjk d~ (j kt A l -* 

A ■ A 

(41) 


^Note that it is possible for A ■ A = 0 without A being zero since the dot product defined here involves 
no complex conjugate. However, this is not so if A is real or purely imaginary. 
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given that cjm is the usual Levi-Civita tensor and the Einstein summation convention on 
repeated indices has been used. 

— I ^ 

We now argue that A remains nonzero when initially nonzero. By taking the gradient of 
(24), it is clear that 

^ = -T + A (42) 

dt v ' 

where the + superscript refers to the transpose of T defined by (39). Therefore, the only 

way for A to be zero in finite time is if the elements of the tensor T blow up in finite time. 

If u s (x, t) is smooth for all times and grows at infinity (for an unbounded domain) no faster 

than linearly in x , then clearly T(xo(t),t) cannot blow up in finite time. Thus, A cannot 

be zero in finite time when initially nonzero. Equations (40) and (42) form a closed system 

of equations for A and C, and thus for Vd(x 0 (t),t) and p(x 0 (t),t). This system of ordinary 

differential equations can be solved numerically once u s and its spatial derivatives are known 

in the complex plane. In the limiting case of a singularity approaching the physical domain, 

further simplifications are possible making it feasible to calculate A and C in an approximate 

manner from knowledge of u s and its gradient in the real domain only. 

The motion of complex singularities of vorticity is significant since its approach towards 
the physical domain can be associated with the appearance of small scale structures at 
places in the real domain that are closest to the complex singularities. Furthermore, it is 
clear that the location of initial singularities in the complex plane is related to the data 
specification in the real domain in an ill-posed way; i.e., one can specify data in the real 
domain with arbitrary but finite precision, yet one cannot nail down the singularity location 
in the complex spatial domain to any precision. Yet, if the complex singularities approach 
the physical domain along some trajectories, one can explain the sensitivity of the later time 
evolution in the physical domain to initial conditions. One can conceivably start with a 
random distribution of singularities in the complex domain and generate an understanding 
of where and when in a physical flow one can expect fine scale structures where viscosity is 
important. The full time-dependent calculation will be the subject of a future investigation. 

3. POSSIBLE IMPLICATIONS FOR THE HYDRODYNAMIC STABILITY 
OF STEADY FLOWS 

In order for a flow to be stable, any disturbance to a basic steady flow cannot grow 
without bound as t — » oo. A flow will be considered absolutely unstable here if there is a 
continual growth of disturbance vorticity at some fixed Eulerian point in the flow field. It 
will be considered unstable in a Langrangian sense if there exists a disturbance that grows 
at some real point that moves with the steady flow. 

If we consider a basic steady flow, and associate the steady flow velocity and vorticity with 
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u s (x) and lo s (x), respectively, as described in the last section, we can consider superposition 
of an initial disturbance so that the initial condition is of the form (8) and (9). Our arguments 
in the last section suggest that any point Xo(t) on the singularity surface d(x 0 (t),t ) = 0 will 
move according to (27). Since u s is now a steady flow, (27) reduces to 

dr 

~~ = *.(2 o (0) (43) 

When a singularity continually moves closer towards the physical domain, one can expect 
smaller and smaller scales to be amplified in the physical domain at some physical point in 
the immediate vicinity of the complex singularity. More precisely, in the limit of a complex 
singularity approaching the real domain, its distance from the real domain will be of the same 
order as the small scales generated in the real domain. As discussed before, singularities of 
the form considered here can at best approach the physical domain continually without 
actually impinging it in finite time. 

The relevant question, as far as hydrodynamic stability in the absolute sense is concerned, 
is whether the disturbance vorticity / ; (d(x, t))p(x, t) grows indefinitely with t at some real 
point x. Since d can never be zero in the real domain and must be bounded below by the 
smallest value of d on the real domain at the initial time (which follows from (24)), it is clear 
that the growth of disturbance vorticity will depend on the growth of p in the real domain. 

We will now show that if there is a trajectory Xq ( t) determined from (43) that asymptotes 
to some real fixed as t — > oo, then it is possible to choose a disturbance such that 
p(xo(t),t) grows exponentially in time. 

Without any loss of generality, we will choose x s = 0. Then, according to the given 
condition, there exists an x 0 (t) such that it approaches the origin as t — ► oo. Near the 
origin, i.e. the stagnation point, it is clear that 


u s (x) ~ Tx 


(44) 


where T is now a constant matrix. Then according to (43) and (44), for t > t 0 (where t 0 
is sufficiently large), xq ( t) is determined approximately by 


dx 0 

dt 


= Tf 0 


(45) 


According to our condition of approach to the origin, T must have at least one eigenvalue a 
with a negative real part. We choose x 0 (fo) to be in the eigenspace corresponding to a. If a 
is real, Xo(to) will be chosen to be completely imaginary. Whatever the case, we have 


x 0 (t) = x 0 (toy {t - to) 


(46) 
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Then, if we choose A(t 0 ) = <f to be in the eigenspace of T + corresponding to the eigenvalue 
cr, according to (42), it follows that 

> 

A{t) = ( e ~ a{t - to) (47) j 

From (46) and (47), it is clear that A ■ x 0 is bounded as t — ► oo. Equation (40) now holds 
where the expression (41) for Sjk is given by 

Sjk = T 3k + e jkl Z, j4 (48) 

Clearly, since ^is a left eigenvector of T (i.e. an eigenvector of T + ), it is clear from (48) that 
it is also a left eigenvector of S. Thus the set of eigenvalues of T and S are identical. Since 
the velocity gradient tensor T is traceless (incompressibility constraint), it is clear that the 
sum of eigenvalues of T, and therefore of S, must be 0. Thus, since there is an eigenvalue a 
with a strictly negative real part, there must exist some other eigenvalue A with a positive 
real part. If ( is a nonzero vector of S in the subspace corresponding to eigenvalue A, it 
follows that if C(to) = £, then 

C(t) = <fe A(t -‘ o) (49) 

Thus, clearly there exist disturbances for which C and therefore p(xo(^),i) grows exponen- j 

tially as t — ► oo. Note that the choice made in (49) is perfectly consistent with the condition i 

C- A = 0 that follows from (37), since f and ( are left and right eigenvectors of S correspond- « 

ing to differing eigenvalues <j and A. Our proof that p(xo(t),t) grows exponentially in t when | 

x 0 (t) approaches x s requires that + 0. In the unusual case when there is no eigenvector | 

of T + with this property, a modified analysis is required by choosing a purely imaginary = 

vector £ in the subspace corresponding to the eigenvalues <r and <j (Note that eigenvalues p 

with a negative real part must be complex in this case). Then, the tensor elements Sjk in 
(48) will not be constants but are periodic in time. In that case, we have been unable to ^ 

show in any elementary way that C(t) must grow with time, though we expect this to be r 

the case. I 

Note, however, that for two-dimensional disturbances superposed on a two-dimensional = 

flow with a stagnation point, C(to ) must necessarily be aligned along the axis, which is 
an eigenvector of T corresponding to eigenvalue 0; therefore, it cannot be chosen to be £. 

Indeed, for two dimensional disturbances, the terms on the right hand side of (38), which 
arise from vortex stretching, are absent and C is a constant. | 

From (24), it is clear that at the stagnation point x s , d(x s J ) = d(x s , 0). Thus if d(x s ,0) r 

is chosen sufficiently small, each of /, ///', f 2 f"/f' 2 appearing in (25) is going to be uni- 
formly small. Thus, it seems plausible that as x 0 (<) approaches x s and p(x 0 (t),t) grows p 
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exponentially in time, p(x 8 ,t) (and hence the disturbance vorticity f'(d)p(x 8 ,t)) will grow. 
However, the mathematical proof of this statement is far from obvious for a general flow, 
where p is determined by the nonlinear equation (25). However, in the Appendix, under 
some additional assumptions which appear reasonable, we are able to show that exponential 
growth of p(xo(t),t) does correspond to a growth of disturbance vorticity up = f f {d)p at x s 
for a steady flow where is a constant and up is small enough to permit linearization. 
We expect that the connection actually transcends these restrictions, though w r e have not 
succeeded in proving this. 

Based on the arguments in this section and the results in the Appendix, we present the 
following conjectures. 

Conjecture 1 

A sufficient condition for the absolute instability of a smooth steady Euler flow in a 
periodic or unbounded domain is that a singularity moving in the complex plane according 
to (27) has the property that as t — ► oo, x 0 {t) -* x s for a real x S} where Ira x o (0) ^ 0. 
This condition is equivalent to the existence of a stagnation point x 5 , where the velocity 
gradient tensor T has at least one eigenvalue with a negative real part. In that case, the 
vorticity at x s will grow in time. 

Remark l 

The above conjecture is only a sufficient condition for instability, though not necessary, 
as shown in a later example. This is because only specific forms of disturbances have been 
chosen- namely, those associated with a specific singularity structure in the complex plane. 
Remark 2 

There is no restriction on the size of the disturbance, though there is a restriction on 
its form. Nonetheless, we do not expect this to be useful in uncovering nonlinear instability 
of flows which are linearly stable since the condition for approach of singularities is equally 
valid for linearized equations. 

Remark 3 

Arbitrary disturbances in the conjecture refer to any arbitrary disturbance in three dimen- 
sions, although the steady flow can be two-dimensional. For two-dimensional disturbances 
superposed on two-dimensional flow, the approach of a complex singularity can be associ- 
ated with the growth of the gradient of vorticity only. This can be seen by just taking the 
gradient of the scalar vorticity equation and evaluating it at the stagnation point, which 
leads to dC/dt = TC, where C = Vu(x s ,t), and T is the two-dimensional velocity gradient 
tensor at the stagnation point. 

Remark 4 

This conjecture can be applied to any steady state flow - even ones that are numerically 
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determined - once the deformation tensor T at a stagnation point is obtained. 

Conjecture 2 

A steady flow is unstable in the Langrangian sense if there exists a choice of initial 
conditions x(0), y( 0), A(0) and C(0) such that the solution to (29), (30), (40) and (42) have 
the feature that as t — > oo, y — » ON and C(t) — > oo. 

In the following subsections, we present examples of the application of Conjecture 1 to 
specific simple Euler flows. Conjecture 2 appears to be harder to apply; at this point, it is 
not clear to us how this can be applied in a fruitful manner. 

3.1 Elliptical Flows 

We first consider the case of plane strain in an infinite flow domain with the velocity field 


u 3 (x, 0) = (axi, — ax 2 ,0) 


(50) 


with a > 0. Then according to (27), a singularity at xq — {xi^x 2} x 3 ) moves according to 


Hence 


dx i 
dt 


= ax\ 


dx 2 

dt 


—ax 2 


dx 3 

dt 


-0 


Xi = Xi(0)e a ‘ , x 2 = x 2 (0)e at , x 3 = x 3 (0) 


(51) 

(52) 

(53) 

(54) 


from which it is clear that if we choose Re Xj(0) = 0 for 1 < j <3 and Im Xi(0) = 

Im ^3(0) = 0, then Xo(t) —> 0 as t 00. According to Conjecture 1, the flow is going 
to be unstable in the absolute sense. This is in accordance with established results from 
hydrodynamic stability theory (see Pearson 12 , Lagnado et al 13 and Bayley 14 ). 

The relation of the approach of a complex singularity to the appearance of increasingly 
small scale structures is highlighted by the following perturbation of our straining flow: 


u(x, 0) = {ax 1 + /(x 2 ), -ax 2 , 0) 


(55) 


for an arbitrary nontrivial differentiable function /. In this case, there is a simple exact 
solution** given by 

u(x,t) = [axi + e~ at f(x 2 e Qt ), — ax 2 , 0] (56) 

II Recall j?o(0 — x{t) + iy(t) 

"""This is probably well known; however, we could not find a specific reference 
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with the corresponding two-dimensional vorticity 

« = -/'(*>£"') ( 57 ) 

Note that if the function / has a singularity in the complex plane, it approaches the physical 
domain as t — * oo. For instance if 

fM = - 7 —=; ( 58 ) 

V 1 + x 2 

then from (55) and (56), it follows that the singularity of u> and u is located at X 2 = i e~ al 
which approaches zero as t — * 00 . In this case, with a two dimensional disturbance super- 
posed on a two dimensional flow, the approach of a complex singularity does not correspond 
to growth of vorticity; however, as expected from Remark 3, the gradient of vorticity grows 
exponentially with time. The solution form (57) in the physical domain clearly corresponds 
to transfer of energy to smaller scales near X 2 = 0 , though the energy contained at any 
specific scale eventually decreases with time. The similarity variable e at X 2 in (56) could have 
been predicted by looking at the solution (54) in the complex plane. Note that I m x 1 and 
Im X 2 need not go to zero for all initial conditions in order for Conjecture 1 to hold. Indeed, 
there exist initial conditions when singularities move further away from the real domain 
which corresponds to a reverse cascade: energy transfer from smaller scales to larger scales. 
This can also be seen in terms of exact solutions to the Euler equation of the form (56) when 
the initial perturbation velocity involves the second component rather than the first. 

Three-dimensional generalizations of the above straining flow can be made for which the 
velocity is given by 

Us = {otx 1 , 0X 2 , 7 x 3 ) (59) 

where a + 0 + 7 = 0 due to (11). The case of 0 = 7 = — corresponds to the axisymmetric 
expansion/contraction which has been studied extensively in the turbulence literature. As 
7 — » 0, the plane strain case considered earlier is recovered. Here, we show how to apply 
Conjecture 1 in its simpler form. The velocity gradient tensor corresponding to (59) is 

fa 0 0 \ 

T = 0/30] (60) 

^ 007 / 

The eigenvalues of (60) are a, 0 and 7 - at least one of which is negative due to the incom- 
pressibility constraint a + 0 + 7 = 0. Hence, according to Conjecture 1 , the flow is unstable 
in accordance with results from hydrodynamic stability (see Pearson 12 and Townsend 15 ). We 
will apply Conjecture 1 to the remaining elliptical flows to be considered with this simple 
eigenvalue test. 
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Plane strain plus rotation in an infinite flow domain corresponds to the base flow 


(cvart — fhrijfhci — cvx2,0) 


This physically corresponds to a plane strain (with strain rate a) combined with a solid body 
rotation (with angular velocity fi) in the plaine of the strain. The velocity gradient tensor 
corresponding to this flow is 

/ a -U 0 \ 

T = fi -a 0 (62) 

V 0 0 0 / 

The eigenvalues of (62) are 0, y/a 2 — ft 2 and — \Ja 2 — Q 2 and, hence, one of them has a 
negative real part if a 2 > fl 2 . Consequently, according to Conjecture 1, this flow will be 
absolutely unstable when 

o' 2 — fl 2 > 0 (63) 

This result is consistent with recent stability analysis of elliptical flows (see Bayley 14 and 
Craik and Criminale 16 ). Since Conjecture 1 only provides a sufficient condition for instability, 
it does not formally apply when fl 2 > a 2 . Nonetheless, it is interesting to note that fl 2 > a 2 
is actually an established condition for the stability of this flow (see Bayley 14 ). 

We now consider the base flow 

u s = (Sx -2 — 2flx2,2flx \ , 0) (64) 

which serves as a model for uniform shear flow subject to a system rotation with angular 
velocity fl (see Batchelor 11 ). Uniform shear flow is recovered in the limit 0 — » 0. The 

velocity gradient tensor T here is given by 

0 5-2fi 0\ 

T= 2U 0 0 (65) 

0 0 0 / 

The eigenvalues of (65) are 0, \j2fl(S — 2fl) and --^2fl{S — 2fl). Hence, according to Con- 
jecture 1 , this flow is unstable if 

2fl(S -2fi) > 0 (66) 

since this condition guarantess the existence of a real negative eigenvalue. In terms of the 
rotational Richardson number Ri = — (l — 2^ first introduced by Bradshaw 17 , it is 
clear that according to Conjecture 1, the flow is unstable if 

Ri < 0 (67) 

This result is identical to that obtained previously by more complicated linear stability 
analyses (see Lezius 18 and Lezius and Johnston 19 ). 
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For a uniform shear flow, Ri = 0 and, therefore, Conjecture 1 is not useful in determining 
the stability or instability of the flow. The same is true for the more general case of plane 
shear flows where 

a. = (/W. o.o) < 68 ) 

Of course, it is well known that uniform shear flow is weakly unstable; the kinetic energy of 
infinitesimal disturbance grow linearly with time 13 ’ 15 . The necessary and sufficient condition 
for a general plane shear flow (68) to be unstable is for /'(■I’l ) to change sign - the classical 
inflexion point theorem of Rayleigh. 

3.2 Kelvin-Stuart Cat’s Eye 

This class of exact steady solutions to the Euler equation has been presented by Stuart 20 . 
Without any loss of generality, his solution for the stream function 'P can be written as 


«P(x,,x 2 ) = ^ In [ cosh (2 x 2 ) + p cos (2 x,)] (69) 

where 0 < p < 1 is a constant characterizing the family of solutions. The velocity field u s 
in this case is given by 

sink 2x 2 P sin 2xi 

cosh 2x2 + P cos2x\ ’ cosh 2x 2 + p cos 2xi 




Within a period 7 r in x u there are two distinct stagnation points at (0, 0) and (|, 0). At 
(|,0) the velocity gradient tensor T is given by 


( 0 

-J£- 

1 -P 

0 


r~ 0\ 

i-p 


0 

0 


0 
0 ) 


( 71 ) 


The eigenvalues of this tensor are 0, 2i/9 1 ^ 2 / ( 1 — p) and —2ip l ^ 2 /(\ — p), none of which have 
a negative real part. Conjecture 1 is therefore not applicable and nothing can be deduced 
about the stability from the existence of this stagnation point. However, at the stagnation 
point (0, 0), the velocity gradient tensor is 


T = 


/ 0 

2 P 

1 +P 


\ 0 


T+p °\ 
0 0 

0 0 / 


(72) 


The eigenvalues of this tensor are 0, — 2p 1/,2 /(l + p) and 2p ly,2 /(l + />), one of which is real 
and negative. Hence, according to Conjecture 1, the Kelvin-Stuart Cat’s Eye is unstable. 
This agrees with the findings of Pierrehumbert and Widnall 21 . 
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4. CONCLUDING REMARKS 


We have considered the motion of a specific type of singularities of the incompressible 
three dimensional Euler equation in the complex spatial domain that correspond to a smooth 
variation in the real domain. Under the assumption that there is a class of initial conditions 
where no new singularities form in the real or complex domain other than what is present 
initially, we have derived an equation (namely, Eq. (27)) for the trajectory of a singularity 
and have shown that such a singularity cannot impinge the physical domain in finite time. 
Nonetheless, we expect that the approach of these singularities towards the physical domain 
will correspond to the formation of small scale structures therein, where the neglected viscous 
effects would be important. The sensitivity of these small scale features to initial conditions 
in the physical domain can then be traced to the problem of determination of the initial 
singularity location in the complex plane from finite precision data in the real domain. 
Furthermore, from (27), we note that the problem of determining the singularity trajectory 
is the complex equivalent of determining fluid particle trajectories in the real domain, once 
an Eulerian flow field is determined. In such cases, “Langrangian chaos” is a well-known 
phenomena. Thus, we speculate that the theory of dynamical systems may be relevant to 
describing chaotic complex singularity trajectories, which could have a bearing on turbulence. 

As far as the hydrodynamic stability of steady flows is concerned, we have introduced 
two new conjectures that connect the complex singularity motion with the growth of a 
disturbance in the real domain. Admittedly, more theoretical work is necessary to put these 
conjectures on a firm foundation. Nonetheless, the application of Conjecture 1 to specific 
examples that include the class of elliptical flows and the Kelvin-Stuart Cat’s Eye give 
results in agreement with previously determined results. Though no new instabilities have 
been discovered so far, the ease and simplicity with which these results have been found is 
certainly encouraging. Indeed, as mentioned before, Conjecture 1 can even be be applied to 
a numerically determined smooth steady Euler flow provided there is at least one stagnation 
point where the velocity gradient tensor has an eigenvalue with a negative real part. 

So far, we have only considered the stability of Euler flows with no boundaries. We do 
not expect boundaries to alter the statement of Conjecture 1 - provided that the stagnation 
point is in the interior of the domain - since the instability mechanism associated with the 
approach of a complex singularity is associated with small scale disturbances in the real 
domain. This issue will be the subject of a future study. 
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APPENDIX 


In Section 3, it has been shown that if the velocity gradient tensor T at a stagnation 
point x s (taken to be zero without any loss of generality) has an eigenvalue a with a negative 
real part, then there exists a singularity trajectory approaching the stagnation point such 
that, for sufficiently large t , the singularity trajectory x 0 (t) asymptotes to 

x Q (t) = x 0 (to)e <r< ‘ t ~ t °' > Ml) 


In that case, it was shown that p(xo(to)^to) can be chosen such that for sufficiently large 
p(xo(t), t) asymptotes to 

p(x o(t), t) = e x ^p(x 0 (t 0 ), to) M2) 


where A is an eigenvalue of T with positive real part. 

The purpose of this Appendix is to show that the exponential growth of p in time (as 
in (A2)) at the moving singular point Xo (t) that approaches the origin (i.e. the stagnation 
point, as in (Al)) corresponds to a growth of the disturbance vorticity up = f{d)p at the 
origin. Here, we only do so for the case when the the velocity gradient tensor at every point is 
T, a constant. We further assume that the disturbance is small enough so that the following 
linearized equations hold 


with 


duo 

~dt 


+ u s ■ Vlo d + u D ■ Vd3 s 


- Co D ■ Vu s + u s ■ Vu D 


Vx«d = LOd 


V ■ u D = 0 


where up = f(d) q. 

We introduce new variables, 


xe-^-'o) = y 
e c(t - to) u D = W 


M3) 

M4) 

M5) 


(A6) 

(A7) 


Then, with the assumption of constancy of Vu s = T (and therefore of tH^), we have 


W t - aW - (tY ■ V Y W + T Y • V y W = TW + u s ■ V Y u D 

Vy X Ud = W 

Vy ■ Up = 0 

where the subscript Y denotes a derivative with respect to the Y components. 


M8) 

m 
M io) 
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Equations (A8) - (A10) are linear equations with no explicit time dependence. Assuming 
that the eigenvectors corresponding to the linear stability operator form a complete set, 
the large t behavior of the solution for a generic initial condition can be expected to be 
dominated by the least stable (or most unstable) eigenvector times e St (multiplied at best 
by a polynomial in t) where 6 is the eigenvalue with the largest real part. The analytical 
continuation of such a solution to the complex plane will have similar time dependence. Thus, 
according to this argument, the time dependence of W at a fixed complex Y = axo(to) can 
only be of the form e Sl (multiplied at best by a polynomial in t). 

However, for the specific initial condition discussed in Section 3, that leads to (A2), it 
follows from the transformations (A6) and (A7) that the time dependence of W(axo(to),t) 
is of the form e (A+ff)i . This growth rate (or decay rate if Re (A + a) < 0)) cannot exceed 
(or have a decay rate less than) e 6i times a polynomial in t. Thus, Re 8 > Re( A + cr) 
implying that for a generic non-zero choice of IT(0, t 0 ), the growth rate of IT(0, t 0 ) will equal 
or exceed e ^ e ( A + ff ) f . From transformation (A7), it follows that u;£)((j, t) will then grow at a 
rate at least equaling e ReA< , implying absolute instability. 
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